Nuclear SchifF moment in nuclei with soft octupole and 

quadrupole vibrations 



N. Auerbach^, V.F. Dmitriev ^, V.V. Flambaum^, 
A. Lisetskiy'^, R.A. Sen'kov^, and V.G. Zelevinsky^ 
^ School of Physics and Astronomy, 
Tel Aviv University, 
Tel Aviv, 69978, Israel 
^ Budker Institute of Nuclear Physics, 
Novosibirsk 630090, Russia 
^School of Physics, 
University of New South Wales, 
Sydney 2052, Australia 
^ GSI, Theory Department, 
64291 Darmstadt, Germany 
^National Superconducting Cyclotron Laboratory and 
Department of Physics and Astronomy, 
Michigan State University, 
East Lansing, MI 48824-1321, USA 



(Dated: February 9, 2008) 



1 



Abstract 

Nuclear forces violating parity and time reversal invariance {V, T-odd) produce V, T-odd nuclear 
moments, for example, the nuclear Schiff moment. In turn, this moment can induce the electric 
dipole moment in the atom. The nuclear Schiff moment is predicted to be enhanced in nuclei 
with static quadrupole and octupole deformation. The analogous suggestion of the enhanced 
contribution to the Schiff moment from the soft collective quadrupole and octupole vibrations 
in spherical nuclei is tested in this article in the framework of the quasiparticle random phase 
approximation with separable quadrupole and octupole forces applied to the odd ^^^~^^^Ra and 
2i7-22ij^j^ isotopes. We confirm the existence of the enhancement effect due to the soft modes. 
However, in the standard approximation the enhancement is strongly reduced by a small weight of 
the corresponding "particle + phonon" component in a complicated wave function of a soft nucleus. 
The perspectives of a better description of the structure of heavy soft nuclei are discussed. 

PACS numbers: PACS: 32.80.Ys,21.10.Ky,24.80.+y 
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I. INTRODUCTION 



The search for interactions violating time reversal (T-) invariance is an important part 
of studies of fundamental symmetries in nature. The manifestations of CP-violation (and 
therefore, through the CPT-theorem, of T-invariance) in systems of neutral K- and B- 
mesons jl| set limits on physical effects beyond the standard model. The main hopes for the 
extraction of nucleon-nucleon and quark-quark interactions violating fundamental symme- 
tries emerge from the experiments with atoms and atomic nuclei, see the recent review P| 
and references therein. For example, the best limits on V, T-odd forces have been obtained 
from the measurements of the atomic electric dipole moment (EDM) in the ^^^Hg P| and 
^^^Xe yi nuclei. As we know from past experience with P-odd forces, see the review article 
P], there are powerful many-body mechanisms in heavy atoms and nuclei which allow one 
to expect a significant amplification of effects generated on the level of elementary interac- 
tions. There are also suggestions for using possible molecular and solid state enhancement 
mechanisms Q, • 

Theoretical calculations of atomic EDM proceed through the nuclear Schiff moment S 
since the nuclear EDM is shielded by atomic electrons J^. The Schiff moment produces the 
V, T-odd electrostatic potential that, in turn, induces the atomic EDM. The expectation 
value of the vector operator S in a stationary nuclear state characterized by certain quantum 
numbers of angular momentum, JM, is possible only for J 7^ owing to the requirements 
of rotational invariance. Since all even-even nuclei have zero ground state spin, we need to 
consider an odd-A nucleus. Furthermore, the non-zero expectation value of a polar vector S 
requires parity non-conservation in a nucleus; in addition, being proportional to the T-odd 
pseudoscalar ((S ■ J)), this expectation value reveals the violation of T-invariance. 

A reliable evaluation of the nuclear Schiff moment should include the estimates of renor- 
malization effects due to "normal" strong interactions inside the nucleus. The core polar- 
ization by the odd nucleon is important, especially in the case of the odd neutron, as ^^^Hg 
and ^^^Xe. Calculations show that the resulting configuration mixing, de- 



pending on details of the method, may change the result of the independent particle model 
by a factor of about 2. A possibility of using accidental proximity of nuclear levels with the 
same spin and opposite parity was pointed out in Refs. In such approaches, possible 

coherent enhancement mechanisms are usually not considered. 
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The statistical many-body enhancement of parity non-conservation in the region of the 
high level density of neutron resonances was predicted theoretically (see e.g. reviews ^, 



15l | and references therein). The existence of such enhancement is now well documented 



experimentally 



lq|. The simultaneous violation of parity and time reversal invariance can 



be enhanced by the firmly established P-violation due to P-odd T-even weak interactions. 
The idea of a possible role of static octupole deformation 

H Q, Q exploited the 

parity doublets which appear in the presence of pear-shaped intrinsic deformation of the 
mean field. The doublet partners have similar structure and relatively close energies so 
that they can be more effectively mixed by P-odd forces. The Schiff moment in the body- 
fixed frame is enhanced bein g in fact proportional to the collective octupole moment. The 



microscopic calculations 







2l\ predict a resulting Schiff moment by two-three orders of 



magnitude greater than in spherical nuclei; this enhancement was confirmed in Refs. 



23|. The uncertainties related to the specific assumptions on P, T-odd forces and different 



approximations for nuclear structure are on the order of a factor 2 for the resulting Schiff 
moment. 

It was suggested in Ref. j3| that soft octupole vibrations observed in some regions of 
the nuclear chart more frequently than static octupole deformation may produce a similar 
enhancement of the Schiff moment. This would make heavy atoms containing nuclei with 
large collective Schiff moments attractive for future experiments in search for V, T-violation; 
experiments of this type are currently under progress or in preparation in several laboratories. 
Recently we performed j3] the estimate of the Schiff moment generated in nuclei with the 
quadrupole deformation and soft octupole mode and showed that the result is nearly the 
same as in the case of the static octupole deformation. 

A related idea is explored in the present paper. It is known that some nuclei are soft 
with respect to both quadrupole and octupole modes, see for example recent predictions for 
radioactive nuclei along the N = Z line j2a| . The light isotopes of Rn and Ra are spherical 
but with a soft quadrupole mode and therefore large amplitude of quadrupole vibrations. 
The spectra of these nuclei display long quasivibrational bands j^] based on the ground 
state and on the octupole phonon, with positive and negative parity, respectively. These 
bands are connected via low-energy electric dipole transitions. This situation seems a-priori 
to be favorable for the enhancement of V, T-odd effects. 

Below we examine the question whether the enhancement indeed exists in some spherical 
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nuclei which have both collective quadrupole and octupole modes. The main weak interac- 
tion mixing can be expected to occur between the levels of the same spin and opposite parity 
that carry a significant admixture of the (particle + 2+ phonon) and (particle + 3~ phonon) 
states with the odd nucleon in the same single-particle orbit. In the odd-neutron nuclei the 
Schiff moment is originated by the proton contribution to the collective phonon. A number 
of such nuclei have an appropriate opposite parity level with the same angular momentum 
close to the ground state. The mixing between the states of above mentioned nature can 



imitate the evo 
with Refs 



Lie CVU 

3. 



ution of the spherical nucleus to the deformed pear-shaped state. In analogy 



21|, this might lead to enhancement of Schiff moments. 



We will show that the effect does exist but, in the framework of the standard quasiparticle 
random phase approximation (QRPA), it is essentially compensated by the strong spreading 
of the single-particle strength that significantly reduces the weight of the desired configura- 
tion in the odd-nucleus ground state. We concentrate here only on physics connected to the 
soft vibrational modes. At the next stage the task of theory should be to go beyond the RPA 
and to combine all effects generated by V, T-violating interactions, including the influence 
of the core polarization and weak interaction admixtures to the phonon structure which 
are not discussed on the present work. Such effects should be included for the quantitative 
answer; however they are not expected to lead to a qualitative enhancement. 



II. COLLECTIVE SCHIFF MOMENT IN SPHERICAL NUCLEI 

Consider an odd-A nucleus with two close levels of the same spin J and opposite parity, 
ground state |g.s.) and excited state \x). The energies of these states are Eg,s. and E^, 
respectively. Let be a {V, T)-odd interaction mixing nuclear states. Assuming that the 
mixing matrix elements of S and W are real, we can write down the Schiff moment emerging 
in the actual mixed ground state as 

g ^ o (g.s.|iy|3:)(x|S|g.s.) 

However, as it was explained in Ref. j^, in the case of mixing of close single-particle 
states one should not expect necessarily a large enhancement. For example, in a simple 
approximate model, where the strong nuclear potential is proportional to nuclear density 
and the spin-orbit interaction is neglected, the matrix element (g.s.|W^|x) contains the single- 
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particle momentum operator. This matrix element is proportional to {Eg,s. — E^), so that the 
small energy denominator cancels out. As mentioned above, the collective Schiff moments 
in nuclei with static octupole deformation may be by 2-3 orders of magnitude stronger than 
single-particle moments in spherical nuclei. 



Consider the 
spherical nuclei 



bl 



owing example of a mechanism generating a collective Schiff moment in 
. Let an odd-A nucleus have the unperturbed ground state of spin J 
built as a zero spin core plus an unpaired nucleon in the spherical mean field orbit with 
angular momentum j = J. The interaction between the odd particle and vibrations of the 
core causes an admixture of a quadrupole phonon to the ground state if the nuclear spin 
J > 1/2: 

|g.s.) = ao|j = J)+a2|[j®2+]j). (2) 

An opposite parity state with the same spin J can be formed by coupling an octupole phonon 
to the ground state particle, 

\x) = \[j®3-]j). (3) 

A (P, T)-odd Schiff moment (P) can mix the states |g.s.) and \x). To first order in the 
non-relativistic nucleon velocity p/m, the ("P, T)-odd interaction can be presented as 



^^^i^iVabCTa - VbaO-b) ■ aS{ra - n) + r]'^f,[cra X CtJ ■ { (p„ - Pfo) , 5(r„ - T;,) } j , (4) 

where { , } is an anticommutator, G is the Fermi constant of the weak interaction, m is the 
nucleon mass, and cra,b, Va^b, and Pa,b are the spins, coordinates, and momenta, respectively, 
of the interacting nucleons a and b. The dimensionless constants rjab and r/^f, characterize 
the strength of the {V, T)-odd nuclear forces; in fact, experiments on measurement of the 
EDMs are aimed at extracting the values of these constants. 

In the case of static deformation, in the "frozen" body-fixed frame the intrinsic collective 
Schiff moment ^intr of the deformed nucleus can exist without any {V, T)-violation. The 



estimate found in Refs. 



2l| gives 



5i„tr ^ — ^ eZR'P^Ps = Ointr/?2 , (5) 

207rv35 5v35 
where P2 and are the static quadrupole and octupole deformation parameters, respectively, 
and Ointr is the static octupole moment. Of course, in the space-fixed laboratory frame, the 
nucleus has definite angular momentum rather than fixed orientation, and this makes the 
expectation value of the Schiff moment to vanish in the case of no ("P, T)-violation. 



The relation (jSJ is expected to hold j25| for the dynamic quadrupole and octupole de- 
formations in systems with spherical equilibrium shape if the effective dynamic deformation 
parameters and (measured by the multipole transition strength for the phonon states 
in the adjacent even-even nucleus) have a magnitude similar to that in deformed nuclei. 
This can be the case under an assumption that the ground state of the considered odd-A 
nucleus can be well approximated by an unpaired particle with an admixture of a single 
phonon of a corresponding multipolarity. It is therefore important to try to perform a de- 
tailed calculation and learn more about the validity of such estimates. Below we present 
such calculations in the framework of the conventional QRPA for quadrupole and octupole 
phonons. 

III. QRPA PHONONS AND PARTICLE-PHONON COUPLING 

We start with the calculation of quadrupole (J'" = 2+) and octupole (J'^ = 3~) phonon 
states in heavy even-even nuclei based on a conventional model Hamiltonian, 



The first term, ff„_p , describes the single-particle mean field. The following Woods-Saxon 
parameterization j28| was used in calculations: 



are the strength parameters for the central and spin-orbital potentials, respectively; R = 
1.3 A^^^ fm and a = 0.7 fm are the radius and diffuseness parameters. The second term of 
eq. (jni), impair, is the pairing interaction with the strength 



H — -f^s-p + -f^pair + -f^phon- 



(6) 




and 




17.9 + 0.176 (N-Z) 
A 



for protons and 



n 



18.95 - 0.078 (N-Z) 
A 
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for neutrons. The BCS formalism was used that yields a quasiparticle basis and correspond- 
ing Bogoliubov transformation coefficients uj and Vj. The last term, 

-f^phon = ^ ^xQif^Qxfj,, (7) 

presents RPA phonons obtained with a separable multipole-multipole interaction. The build- 
ing blocks of the model are the two-quasiparticle RPA phonons, 

Ql = lll K.[4«]Jv - (-)'-^5LK«..]a-m] , (8) 
where Aj-^^j^ and -Bjijj are the forward and backward phonon amplitudes, 

where Z{X) is a normalization factor, f'^{jij2) =< Ji|k'^^A||j2 >, and the coherence factors 
of the Bogoliubov canonical transformation are ^j^^ = U1V2 ± f 1M2 and r^J^-* = U1U2 ± ViV2- 
The phonon frequencies ux are the roots of the characteristic RPA equation for each 
multipolarity A, 

where xW is the strength parameter and e(ji,j2) = ^(ji) + ^(j2) is the unperturbed two- 
quasiparticle energy. Solving these equations one obtains the energies of the phonons and 
internal structure of the phonon operator (jH)) hidden in the amplitudes of different two- 
quasiparticle components. The values of the strength parameters, x(2) and x(3), were 
chosen to reproduce the excitation energies of the 2^ and 3i states in the neighboring 
even-even nuclei. 

We are interested in the evaluation of the Schiff moment for the odd-neutron nuclei. 
Therefore the second step of calculations is the solution for the odd-A nucleus with the 
even-even core excitations described above. If we neglect the (quasiparticle+two-phonon) 
components in the wave function of an excited state of an odd-A nucleus, then the corre- 
sponding wave functions have the following form: 

^nir) = Cjn(a^j,,6jj + Y,DjFH ® Q+(A)]jAf)^o, (11) 
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TABLE I: Phonon energies in even-even nuclei and the corresponding parameters of the multipole- 
multipole interaction. 



Nucleus 


u;2, MeV 


u;3, MeV 


X(2), 10-4 


X(3), 10-5 


|g^Rai28 


0.69 


1.40 


8.1 


1.2 


|g^Rni3o 


0.46 


1.11 


6.9 


1.2 


Ig^Raiao 


0.39 


0.79 


6.9 


1.2 




0.32 


0.84 


5.9 


1.1 


gg°Rai32 


0.18 


0.47 


5.9 


1.1 


gg°Rni34 


0.24 


0.66 


5.3 


1.1 



where the operator creates a quasiparticle with respect to the ground state wave function 
\l/o of the even-even nucleus. The energy -E„(J) of the rfi^ state with angular momentum J 
in the odd-mass nucleus, the amplitudes of the quasiparticle-phonon components, 



y i2J+l)2ZiX) Eij)+Lu,-E^' ^ ' 



and the amphtude of the single-quasiparticle component, 

are obtained from the solution of the secular equation [s^ : 

e{J) -E^ = J2 {D'y-y m +UJ,- E^] . (14) 

We performed the calculations for several isotopes of radium and radon. The lighter even- 
even isotopes, A = 216 — 218, are considered to be still spherical; the ones with A = 220 
are transitional and the heavier isotopes A = 222 — 224 acquire more pronounced deformed 
structure. In the six even-even nuclei treated here, see Table 1, the 2+ and 3^ states were 
experimentally observed with excitation energies uj2 and listed in Table 1. The values of 
the coupling constants x(2) and x(3), also shown in this Table, were found by fitting the 
lowest RPA 2+ and 3^ states to the corresponding experimental excitation energies. 

The single-particle states were calculated in a spherical Woods-Saxon potential with the 
parameters indicated above. For heavier isotopes this might not be an appropriate approx- 
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imation; the purpose of our calculations for these nuclei was to examine the systematics of 
the results in the cases when the phonons are lowered in energy. 

In Table II we show the composition of the ground state wave functions and their suitable 
parity partners in the six odd-neutron nuclei considered in this work. One can see from 
the wave functions that while the dominant component for negative parity is the 3^ ® (79/2 
configuration, the corresponding component 2"'' $5(79/2 is of the order of 10-30%. The spacings 
between the ground states and their parity partners turn out to be of order 1-2 MeV, Table 
I. The mixing of these two components will contribute the most to the Schiff moment. We 
should note here that the calculated spectra of the odd-neutron nuclei do not agree well 
with experiment even in the cases of A = 217, 219. This indicates that the model used for 
the description of coupling of the odd particle with the core vibrations is too restrictive. 



IV. THE SCHIFF MOMENT 



Having determined the wave functions one may proceed with the calculation of the Schiff 
moment according to Eq. (^, where we will use the notations \Jg) and \Jg'^) for the 
ground state of the odd-neutron nucleus and for the first excited state of opposite parity, 
respectively. 

Here W is the {V, T)-violating nucleon-nucleon interaction given by Eq. and 
{Jg'^Jg\S\JgJg) is the matrix element of the Schiff operator S, 



3 ^2 

'i 3 ch'* 



(16) 



for the maximum projection M = Jg of angular momentum, and r"^^ is the mean square 
charge radius. 

Since only the proton components of the phonons contribute to the matrix element of 
the Schiff moment, this matrix element for the case of the odd neutron can be written as 

{jfM = Jg\S.\jjM=Jg) 
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TABLE II: Structure of the ground ( J'^) and excited state with the same spin and opposite parity 
(J""^) in odd-neutron nuclei (probabilities of various components in the wave functions). 



88 "-^129 


9/2+ 


259/2 0.71 
2g9/2®2+ 0.14 
1J15/2 <^ 3- 0.08 

3^5/2 + 2+ 0.05 


9/2- 


1/19/2 10-4 
2^9/2 (8) 3- 0.999 




9/2+ 


2^9/2 0.73 
2^9/2 2+ 0.11 
lji5/2 ® 3- 0.07 
34/2(8)2+ 0.07 


9/2- 


lhg/2 2 • 10-4 

2C/9/2 <8) 3- 0.999 


88 i^aisi 


7/2+ 


2ff7/2 0.43 
lin/2 <8 2+ 0.27 
2c/7/2®2+ 0.11 
299/2(8)2+ 0.08 
3^3/2(8 2+ 0.09 


7/2- 


2/7/2 0.02 
2^9/2 <8) 3- 0.98 


86 Rni33 


5/2+ 


3^5/2 0.48 
2ff9/2<8 2+ 0.27 
34/2 (8 2+ 0.12 
4si/2(8)2+ 0.09 


5/2- 


2/7/2 0.01 
2C/9/2 (8) 3- 0.98 


gg-^Raiaa 


5/2+ 


34/2 0.48 
2^9/2(8 2+ 0.25 
34/2 (8 2+ 0.13 
4si/2«)2+ 0.10 
34/2(8)2+ 0.03 


5/2- 


2/5/2 0.04 
2^9/2 (8 3- 0.86 
liii/2 (8) 3- 0.09 


221 R„ 
86 tini35 


7/2+ 


2^7/2 0.44 
liii/2 (8 2+ 0.26 
2c/7/2 (8 2+ 0.12 
34/2(8)2+ 0.12 
2^9/2(8 2+ 0.04 


7/2- 


2/7/2 0.02 
2C/9/2 8) 3- 0.98 
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TABLE III: Mixing matrix elements of weak interaction (m.e. W), of the Schiff moment (m.e. S), 
and final value for the Schiff moment, S, in the ground state of the odd-neutron nucleus. 



Nucleus 


E+, MeV 


MeV 


m.e. W, 7] ■ 10 ^ eV 


m.e. S, e-fm^ 


S, T] ■ IQ-^e-fm^ 


gg'^Rai29 


0.09 


2.44 


0.3 


-0.1 


-0.03 


ie^Rnisi 


-0.12 


1.80 


0.1 


-0.1 


-0.01 


gg^Raisi 


0.49 


1.47 


-1.3 


-0.1 


0.30 




-0.09 


1.42 


0.2 


-0.1 


-0.03 


gg^Rai33 


-0.77 


1.05 


0.2 


-0.2 


-0.07 


||^Rni35 


0.24 


1.16 


-0.5 


-0.1 


0.06 



where D = C ■ D - is the amplitude of the (quasiparticle plus phonon) component, and the 
sum over A and A' is in reality reduced to one term with A'^ = 2"*" and A''^ = 3". The matrix 
element of the Schiff operator between the quadrupole and octupole phonon states in the 
even nucleus has the following form in terms of the RPA amplitudes (0): 

(2^115113-) = V35Y: il0ill^llj2) I ' ' ' I (A^A^^ + bSSbS^)- (18) 

ilj2J3 [ Jl J2 J3 J 

The results of the calculations are shown in Table III. 



V. MATRIX ELEMENT OF WEAK INTERACTION 

The calculation of the weak interaction matrix element appearing in Eqs. (Q) and ()15p is 
elaborate and contains a number of contributions. We collect these calculations more in de- 
tail in Appendix. Here we present only the final expressions for the important contributions 
and numerical results. Table III. 

In this work we account only for the direct contribution of the weak interaction, Eq. (^, 
to the nuclear mean field for a particle b = p,n at a point r = rn, 

where v'^ is the occupancy of the single-particle state c in the core. We expect this term 
to be the largest since all particles in the core add coherently in Eq. |T^. The weak 
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matrix elements are typically of the order 10~^ rj eV. The main contribution comes from 
the mixing between the (quasiparticle plus phonon) component of the excited state and the 
pure quasiparticle component of the ground state. 

The reduced matrix elements (2+115113^) in even nuclei are of the order (l-2)e ■ fm^, and 
the matrix elements of the Schiff operator between the ground state and its parity partner 
turn out to be around 0.1 — 0.2 in units of e ■ fm^. The results for the Schiff moment in 
this model are within the limits (0.04-0.6)7] ■ 10~®e ■ fm^, i.e. of the same order as in pure 
single-particle models in even nuclei. The single-particle contribution unrelated to the soft 
III, 12, 3 of the same order of magnitude and it should be added. 



modes 



VI. DISCUSSION 

How can one reconcile this result with the large enhancement found in the calculations 
for the nuclei with static quadrupole-|-octupole deformation? We need to stress that in those 
calculations the deformations are introduced explicitly from the beginning (in some calcula- 
tions, only the quadrupole deformation is taken at a starting point while the octupole part 
is introduced via coupling of a particle to the 3~ phonon of the deformed core). Apparently, 
it is essential to start from a particle moving in a quadrupole deformed field in order to find 
large enhancement factors for the Schiff moment. The attempt to reproduce this effect by 
having in the ground state of the odd-A nucleus an admixture of a particle plus a single 2"*" 
phonon can succeed only in limiting cases. 

In the case of static quadrupole, P2, and octupole, /?3, deformations, the weak matrix 



element is oc /^s, the Schiff moment operator oc I32(3s/ AE± as discussed in j20|, |2l|. Although 
the small denominator AE± ~ 50 keV was the main source of enhancement, it was also 
important to have noticeable values of the deformation parameters. Since in dynamical 
models with low phonon frequency u the dynamic deformation parameters (3 oc 1/uj, one 
may expect a considerable enhancement a.t uj 0, close to the phase transitions to static 
deformation. Such hmiting cases will be discussed below. 

The effects of quadrupole deformation can be mimicked by a high degree of coherence in 
a quadrupole mode. Similarly to that, the transition of a spherical nucleus to an octupole 



deformed state was discussed in the RPA framework in the past 



It was found that even 



for realistic interactions the energies of the lowest 3 RPA root collapse in certain areas 
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of the periodic table. The stabihty of the nucleus against variation of a certain collective 
variable Q can be expressed in terms of the static polarizability defined for instance through 
the dependence of constrained Hartree-Fock energy E/Ao\i the constraining term aQ, 



2 



d^{E/A) 



J a=o 



(20) 



In the RPA, this quantity can be calculated as 



P = 2 — i (21) 



by relating it to the m„i moment (the inverse energy weighted sum rule), 



^ I^^HVIWI ^ (22) 



with LOn being the RPA energies corresponding to excitation from the ground state |0) to 
the excited state |n) with quantum numbers dictated by the symmetry of the operator Q. 
As polarizability increases the system becomes less stable against deformation characterized 
by the collective variable Q] a negative value of P indicates a situation in which a collapse 
has occurred (in the RPA framework). 

In order to take advantage of such properties and examine the behavior of the Schiff 
moment in situations close to the RPA collapse that signals a transition from a spherical to 
a pear-shaped nucleus (quadrupole + octupole deformation), we can consider the limit of 
the RPA frequencies and cus going to zero. Introducing formally a small scaling factor y 

t^2,3 Z/t^2,3, y < 1, (23) 

we separate the singular in this limit part of the secular equation (fTUj) as 

The normalization factor Z[X) goes to zero linearly in and the phonon amplitudes 
grow oc As a result of this overcritical coupling, the energy of the ground state in 

the odd nucleus collapses as well (goes to — oo as 1/y^). The wave function components 
flll|l have a constant limit with the single-particle amplitude — 1/2; the same effect was 
observed in for the strong coupling of a quasiparticle with phonons. 

In Table IV the results of the calculation for ^^^'^^^Ra are shown as a function of y- 
scaling. As the collective frequencies go down, the reduced matrix element (2| 15*113) in 
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TABLE IV: Scaling of RPA solutions with the low-lying frequency of collective oscillations in the 
even core: normalization factors Z{X); the reduced matrix element (2||5'||3) of the Schiff moment in 
the even nucleus; energies and of the ground and excited state in the odd nucleus, in MeV; 
the matrix element (m.e. W) of the weak interaction between these states; the matrix element 
(m.e. S) of the Schiff moment in the odd nucleus; and the final result for the expectation value of 
the Schiff moment; the matrix elements of W and S are given in the same units as in Table III. 



Nucleus 


y 


Z{2) 


Z(3) 


(2||5||3) 


E+ 




m.e. W 


m.e. S 


S 


2i9Ra 


1.0 


479 


21696 


1.7 


0.5 


1.5 


-1.3 


-0.1 


0.3 




0.1 


43 


1825 


20 


-3.8 


-1.9 


1.1 


-0.2 


-0.2 




0.01 


4 


182 


195 


-17 


-19 


53 


-0.2 


6.2 




1.0 


324 


18100 


2.2 


-0.8 


1.0 


0.2 


-0.2 


-0.1 




0.1 


31 


1650 


23 


-6.1 


-2.9 


-19 


-0.5 


6 




0.01 


3 


165 


235 


-23 


-21 


-253 


-2.7 


560 



the even nucleus, the matrix element of the weak interaction in the odd nucleus and the 
resulting Schiff moment become large. This is especially pronounced in ^^^Ra, where the 
Schiff moment is strongly enhanced at the collapse threshold. 

In the RPA framework the matrix element of the Schiff moment between the states of 
opposite parity in the odd nucleus can be presented in terms of the amplitudes Djj,, eq. 
(HH), as 

{+\S\-} oc J2 D^^>^^^D^f-\X,\\S\\X,). (25) 

AiAai' 

In the case of ^^^Ra, at the experimental values of the phonon frequencies, y = 1, there 
exists only one large combination of amplitudes, 

\Dl^+\ , |2 = 0.08, \Dl^-\ , ? = 0.98. (26) 

I 99/2:99/2 1 ' I 99/2:99/2 1 ^ / 

Here the excited negative parity state is an almost pure combination 2gQ/2 ® 3~. With 
reduction of frequencies, the amplitudes rapidly decrease going to limiting values of 
|^2(+)|2 ^ g Q]^ g^^^ |£)3(-)|2 ^ 0.02. The wave function (quasiparticle + phonon) turns 
out to be spread over many states, and the enhancement of the Schiff moment due to the 
growth of the matrix element (2 1 15*113) starts with further diminishing frequencies only after 
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the spreading saturates in a given single-particle space. The natural conclusion is that the 
enhancement of the Schiff moment in the presence of soft quadrupole and octupole modes 
is a real physical effect which is however hindered by the complexity of the mixed states in 
the limit of strong coupling between quasiparticles and soft mode phonons. 

However, the ansatz for the wave function of the odd nucleus becomes invalid in 
this limit, and many-phonon components carry a large fraction of the total normalization. 
Leaving the full solution of this problem for the future, we just show that the enhancement 
occurs beyond the RPA framework. 

VII. SIMPLE MODEL BEYOND RPA 

The strong quasiparticle-phonon coupling leads to stationary states which are close to 
the coherent states, or condensates, of the phonon field . This limit is different from that 
used in (jllj) . where only a one-phonon admixture was taken into account. 

We consider the Hamiltonian (for simplicity with one single-particle j-level but the ex- 
tension is straightforward): 



where x\ are the parameters of the particle-phonon coupling, Qx^ phonon operators, and 
Pm'm = ciln^^m' IS the Operator of the single-particle density matrix that satisfies the operator 
equation of motion: 




+ 




(27) 




(28) 



J2 



fi — m m 



^ 3 J 



) 




fm"{Q\, + {-Y-'Qx-,) 



The phonon operators satisfy 




(29) 
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For monopole phonons, A = 0, the exact solution j33| leads to the coherent cloud of 
phonons with the Poisson distribution of the phonon numbers. For A 7^ the analytical 



3J] that allows to separate 



solution does not exist; however, we can use the method of Ref. 
the main effects in the adiabatic limit of slow collective motion. We assume that the even 
core is spherical with low-frequency collective modes, uj^ and ^^3. Then in the odd nucleus 
we have a phonon condensate (Qa^i) that can be found from (j^ . 

This condensate is created by the average odd-particle density matrix {pm'm) that determines 
self-consistently a direction of the intrinsic quantization axis. This direction is arbitrary 
in agreement with full rotational invariance, and we can expand this density matrix over 
multipoles, 

(P^^'^) = $^(-)^+^-™ I ^ ^ ^ (31) 
LA yA -m m' J 

where the spherical function shows the "rotational behavior" of the condensate in the labo- 
ratory frame of the odd nucleus. Then we find the relation between the condensate and the 
density matrix multipoles through effective deformation parameters Px, 

(Qa,) = /SaF,;, /5a = ^^3^ pa. (32) 

As mentioned earlier, the effective deformation parameters induced in the odd nucleus by 
soft modes of the even core are inversely proportional to the corresponding frequencies. (In 
our notations we have Y^^ because this function lowers the projection by fi, in accordance 
with the definition of Qx^.) 

The relation between the states of the odd and even nuclei is established through the equa- 
tion of motion for single-particle operators linearized after the substitution of the phonon 
dynamics by the condensate, 

[am,H]=eam + 2j2^A^ ^ ^ (^m'PxY;:^. (33) 

A/.;m' -m m' j 

In order to diagonalize this set of equations, we are looking for the particle operator in the 
form 

<^^ = Y.^'mkCk. (34) 
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where Ck are particle annihilation operators with projection k onto the axis defined by the 
phonons, and Dl^^. are Wigner rotational functions. 

After simple algebra, this procedure results in the occupation numbers of the particle 
in the self-consistent phonon field, 

{c\cki) = UkSkk', (35) 
and the "Nilsson-type" energies of the odd particle in the body-fixed frame, 

S, = e + y^ ^ ^' ^ ] (36) 

similar to how it would appear as a result of a self-consistent mutlipole-multipole interaction. 
The original lab-frame density matrix (|HT|) is diagonalized in the multipole form, 



47r I J J ^ 

Pa 



■ 2A + 1 V \ -A; A; ' 

In particular. 



where we get the total number of unpaired particles (equal to 1 in our simplified case). 
The model can be advanced further taking into account the non-adiabatic corrections to the 
condensate, analogously to the frequency expansion in (j21j). 

Now we can show that there is a strong enhancement of the Schiff matrix element com- 
pared to the value 5*° = (2 1 15*1113) in the even nucleus. The collective contribution to this 
operator is given by 

Siu = s°Y^i-r^' I ^ ^ ^ ](Ql'Q2, + i-r^''Ql,Q3,')- (39) 



Using the deformation parameters ()32|) for the odd nucleus, we obtain 



^ = 2f3,f3sJ2i-r'' i ^ ^ ^ ] Ys,'Y,_, = - ^ (3,f3,Y,l. (40) 
^ ^ \ -u u -u' / 



-V jJ — /i 

With the self-consistent values of deformation parameters fl32|l 

^i^ 1 0:2X3 



5*° 35^/n UJ2i03 
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P2P3Y,l. (41) 



Consequently, in the odd nucleus we reduced the problem to that for intrinsic deformation 
and corresponding intrinsic Schiff moment, compare eq. (0). In contrast to the RPA, 
the wave function of the odd particle is not fragmented being a self-consistent mean field 
solution in the presence of the condensate. The enhancement due to small frequencies is 
clearly present in (j4H) with no reduction that appeared in the RPA calculation. 

VIII. CONCLUSION 

We presented the calculations of the Schiff moment in odd-neutron nuclei with soft 
quadrupole and octupole vibrational modes present in the even-even spherical core. The 
conventional approach used above is based on the QRPA which allows one to define mi- 
croscopically the structure of collective modes and coupling of the odd particle to phonons 
of those modes. We confirmed the effect of enhancement of the Schiff moment for small 
quadrupole and octupole frequencies. The even-even system close to the onset of deforma- 
tion acts essentially similarly to the statically deformed system, where the enhancement was 
estabhshed earher. 

However, the QRPA supplemented by the standard ansatz for the (particle -|- one phonon) 
states in the odd system, turns out to be inadequate. The enhancement requires strong 
particle-phonon coupling. In this case, the single-particle strength is widely fragmented, 
and a small weight of mixing configurations (particle + quadrupole phonon) and (particle 
+ octupole phonon) with the same spherical orbit for the particle practically cancels the 
enhancement of the Schiff matrix element between 2+ and 3~ phonons. The outcome for 
the Schiff moment in this type of calculation is on the same level of magnitude as for purely 
single-particle states. The QRPA calculations do not predict small energy denominators A± 
which would provide a strong enhancement of the effect. In general, the spectra of odd- 
neutron nuclei in the case of soft transitional even core cannot be reliably reproduced by the 
standard version of the QRPA. The experimental (poorly known) spectra of ^^^'^^^Ra show 
the states, sometimes with unknown spin, of parity opposite to that of the ground state at 
AE < 100 keV. 

We demonstrated the presence of the effect by artificially diminishing vibrational fre- 
quency to the limit when the single-particle spread is saturated. We also sketched the more 
general formalism beyond the RPA that is capable of describing the strong coupling limit 
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for a particle in the field of soft phonons. The self-consistently emerging coherent phonon 
state acts analogously to static deformation, in agreement with the original idea. We hope 
to present corresponding advanced calculations elsewhere. 

We have considered only the lowest symmetric quadrupole phonon in the present paper. 
In the presence of few active protons and neutrons, there is a possibility of low-lying "mixed- 
symmetry" quadrupole phonon in vibrational nuclei that is known to be connected to the 
octupole 3~ phonon by a considerably stronger El strength than the symmetric lowest 
quadrupole phonon |35|]. Therefore one may expect additional enhancement if the coupling 
to the "mixed-symmetry" phonon is taken into consideration. 

The conventional core polarization effects are also important for getting a certain quan- 
titative prediction, although they are not expected to lead to a significant amplification of 
the effect. In the future the full calculations taking into account on equal footing collective 
and single-particle effects, including the P, T- violating corrections to the vibrational modes 
and core polarization, are to be carried out. It is desirable to extend the studies to lighter 
nuclei A ^ 100 with pronounced octupole phonon states and more experimental data 
available. In particular, there are indications for the intrinsic coupling between quadrupole 
and octupole degrees of freedom, see global systematics 13] and recent experimental data 
37[. This may help to understand better the possibility for the enhancement of the Schiff 
moment and to explore physical correlations between the calculated values of the Schiff 
moment and complex nuclear structure. 
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Appendix 

Our wave function of the odd nucleus includes two types of states, pure quasiparticle (QP) 
and quasiparticle+phonon. Then, any matrix element is a sum of four terms. One term 
mixes pure QP components of the wave functions, the other two mix QP and QP+phonon 
components, and the last one mixes two QP+phonon components. It is convenient to 
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w I W I w 

X X X 

(a) (b) (c) 

FIG. 1: Contributions to the weak matrix element of (QP| VF^ean I QP) type- 
illustrate graphically different contributions to the weak matrix element. 
1. QP and QP Mixing 

Here we calculate the weak matrix element between pure QP components in the mean 
field approximation (fT^ : 

MEl = {qP\W„\QP) = (^o|«Jo«i«.«jjV^o)r/i+^W^. = 

Using Eq. (jH)) we can express the products of aj^a^, via phonon creation and annihilation 
operators Q\ and Qx: 

The first term corresponds to Fig. while the second term corresponds to Fig. CJj. The 
sum over A includes all RPA solutions of Eq. (|Tn|l . It is interesting to recalculate this matrix 
element in a different way: 



{ipo\ajo%aiyaj>^\'ipo) = (^o|ai/ajo«]iajJ^o)- 
In this case the matrix element will be expressed as 

MEl = y (_i)^o+J.-Ja^a .a 2Ja^ (+)— 



that corresponds to Fig. Both calculations should give the same result due to the 

completeness of RPA solutions. If we leave only the most collective phonon in the sum, the 
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!W 

X 



(c) 




Jq 



FIG. 2: Contributions to the weak matrix element of (QP| WmeanlQP + phonon) and 
(QP + phonon|T4^inean|QP> types. 



result will be different. In this case we prefer the first way of calculation where the pure 
QP contribution appears explicitly. In fact, this ambiguity is of minor importance since the 
weight of the QP component in the wave function of the excited state is rather small. The 
final result is: 

A 

2. QP and QP + Phonon Mixing 

Now we calculate the weak matrix element between pure QP and QP+phonon compo- 
nents: 

ME2 = (QP|iy^ean|QP + phoUOu) = 

After simple transformations ME2 can be presented as follows: 



Here the first term corresponds to Fig. and the second term corresponds to Fig. |2b- 
Note the different Bogolubov factors that appear here compared to those in the previous 
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FIG. 3: Contributions to the weak matrix element of (QP + phonon| WmoanlQP + phonon) type, 
matrix element. 

3. QP+Phonon and QP Mixing 

Here we calculate the similar weak matrix element between QP+phonon and QP com- 
ponents: 

ME3 = (QP + phonon|#^ean|QP) = 
The resulting expression is: 

Again, the first term corresponds to Fig. and the second term corresponds to Fig. |211. 
4-- QP+Phonon and QP+Phonon Mixing 

There are several terms in this type of the weak matrix element. Some of them are 
presented in Fig. El there are also additional contributions with Ai <-> A2 and Ai = A2. 
Corresponding diagrams for Ai = A2 are similar to those in Fig. [T^,b. 

MEA = (QP + phonon|#mean|QP + phonon) = 
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It is convenient to split this matrix element into several terms, 

{il)o\Qx^aj^ala^al^Q{^\i)Q) = 6j,f,6j,^{ilJo\Qx,QiJipo) + 6j,j^{ilJo\Qx,ala^Qlj'^o)- 

-Sji,i{MQxi(^\2^J.Qx2\^o) - h2u{'ipo\Qx,alajiQlj^o) + {^olQxialal^a^aj^Qljtfjo) . 

The second term in this equation corresponds to Fig. IHll. Here, the weak interaction mixes 
2"*" and 2" or 3~ and S"*" phonon states. Since we do not include 2~ and S"*" phonons in our 
consideration, we omit this term. 
Part 1 : 

Sj^f.Sj2u{i'o\QxiQ{ji'o) SxiX2Vj^j2Wj,j2- 

The result here is given by 

{+) 



ME4.1(JiAi, J2A2) = 6x,x2V:J:J2^J^J2■ 

The corresponding graph is similar to that of Fig. where one should add a parallel 
phonon line and set Ji, J2 instead of Jq, Jq for the QP lines. 
Part 2 : 

-hit^i'ipolQx.al^aj.Qlj'^o) ~ -5j,f,{ilJo\QxAal2^J,,Qx2]\'^o)- 
After simple calculation one obtains 

ME4.2(JiAi, J2A2) = V(2Ai + 1)(2A2 + 1) J^^fi^AJ. 

JaJb 

The second term here corresponds to Fig. and the first term corresponds to Fig. Et- with 
Ai ^ A2. 
Part 3 : 

-Sj2u{MQxi(^l(^JiQx2\^o) ~ -5j2i.(^o|QAi[aiaji,Qlj|V'o)- 

The result is 

.(+) 



ME4.3{J,X,, J2A2) = v/(2Ai + l)(2A2 + l) J^ryitl^^. 



J2 
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[ A2 Jo Ji J ^Jo + i-^ 

The second term here corresponds to Fig. ^jp, and the first term corresponds to Fig. |2b 
with Ai ^ These terms are similar to those calculated in Part 2. 
Part 4 : This is the most complicated part of the calculation. 



Here 

^ _ ^_-|^^A3+/i3+A4+/i4 

(V'oIQasQaiQI^QIJV^o) = 5aiA2^A3A4 + '^AiA4'^A2A3, 

and we use the RPA approximations 

After some calculations we obtain 

Mi?4.4(JiAi, J2A2) = v/(2Ai + l)(2A2 + l)5^r/i+i^T^j„j, 

ah 

Aj//i 

The second term corresponds to Fig. Eb, and the first term corresponds to Fig. Efc with 
Ai ^ A2. The last term corresponds to Fig. ^ with the additional parallel phonon line and 
Ji, J 2 instead of Jo, Jg for the QP lines. 
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